High-energy neutrinos from reverse shocks in choked and successful relativistic jets 
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Highly relativistic jets are a key element of current gamma-ray burst models, where the jet 
kinetic energy is converted to radiation energy at optically thin shocks. High-energy neutrinos are 
also expected, from interactions of protons accelerated in the same shocks. Here we revisit the 
early evolution of a relativistic jet, while the jet is still inside the star, and investigate its neutrino 
emission. In particular we study propagation of mildly relativistic and ultrarelativistic jets through 
a type lb progenitor, and follow reverse shocks as the jets cross the star. We show that protons 
can be accelerated to lO^-lO^ GeV at reverse shocks, and efficiently produce mesons. The mesons 
experience significant cooling, suppressing subsequent neutrino emission. We show, however, that 
the neutrino yield from the reverse shock is still reasonably large, especially for low-luminosity and 
long-duration jets, where meson cooling is less severe. We discuss implications of our results in the 
context of neutrinos from choked jets, which are completely shock heated and do not break out of 
the star. From a choked jet with isotropic equivalent energy of 10^^ erg at 10 Mpc, we expect ~20 
neutrino events at IceCube. 

PACS numbers: 97.60.Bw; 98.70.Rz; 95.85.Ry. 
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I. INTRODUCTION 

The leading model of gamma-ray bursts (GRBs) in- 
volves a relativistic (Tj ^ 100, where Tj is the jet Lorentz 
factor) fireball jet, where the observed gamma rays are 
produced by radiation from Fermi-accelerated electrons 
in optically thin shocks (for reviews see e.g. [l|, H, [3])- 
In the most popular model for the more common long- 
duration GRBs, the so-called coUapsar model, the core 
of a massive star collapses to a black hole or neutron 
star, driving a highly relativistic jet which breaks out of 
the star Within this scenario, high-energy neutri- 

nos from photomeson interactions of accelerated protons 
have been studied [1, 0, i, i, [M [H, [11 0, El 111, [S- 
With 1 km'^ size neutrino detectors such as IceCube |17| 
and KM3Net [l^ being constructed, we are entering an 
age when we can test various predictions for high-energy 
neutrinos. Waxman and Bahcall predicted Ci, > 10^ TeV 
neutrinos from external reverse shocks |9|, and similar 
predictions have been made in Refs. [ll.]l4l| for the for- 
ward shock. In addition, internal shocks can occur while 
the relativistic jet is still in the star. Meszaros and Wax- 
man predicted a 5 TeV neutrino precursor burst 
from such a set up [161]. It is expected that a precursor 
neutrino burst is present even if the jet does not success- 
fully break out of the star, producing a choked jet that 
is dark in gamma rays and bright in neutrinos. 

Inferring from the observed rate of GRBs, highly rel- 
ativistic {Tj > 100) jets occur in perhaps < lO^'^ of 
core-collapse supernovae fl^. However, the jet signature 
may be more common, and a significantly higher portion 
of supernovae may be endowed with slower, mildly rel- 
ativistic jets; according to late-time radio observations, 
perhaps as high as a few % [19, .20, 21, 22, 23]. The re- 
cent detection of low-luminosity GRB 060218 also sug- 
gests from detection rates that low-luminosity GRBs, 



with mildly relativistic ejecta, are more common than 
conventional high- luminosity GRBs by factors as large 
as 10^ [H, m 111 m m, [ia, US. The occurrence rate 
is high enough that despite their lower luminosity, their 
contribution to the neutrino background is comparable 
to or larger than conventional GRBs [11 [sj] • 

In addition to being more frequent, mildly relativis- 
tic jets are more baryon rich, which work positively for 
neutrino emission. Indeed, the detectability of neutrinos 
from internal shocks of mildly relativistic jets has been 
studied by Razzaque, Meszaros, and Waxman [33] and 
extended by Ando and Beacom [33| (see also [34 \3^). 
In these scenarios, internal shocks occur at rs ~ 10^°- 
10^^ cm, smaller than the typical radius of type lb pro- 
genitors. Accelerated protons interact via proton-photon 
(pj) and proton-proton (pp) interactions to produce 
mesons, which consequently produce neutrinos on decay. 

In this paper we investigate neutrinos arising from pro- 
tons that are accelerated in the reverse shock, while the 
relativistic jet is still propagating inside the star. We 
adopt a type lb progenitor as our model star, but we dis- 
cuss how our results apply to type Ic and II progenitors 
also. The reverse shock is strongly motivated because 
it is a generic prediction of the initial stages of any jet 
moving through the star, regardless of its Lorentz factor, 
and whether it breaks out of the star or not. Therefore, 
we consider two types of supernovae: one is the perhaps 
more numerous supernova containing a mildly relativistic 
{Tj 10) jet, and the other is a supernova containing an 
ultrarelativistic {Tj ^ 100) jet that can ultimately cause 
the observed GRBs. In both cases we consider choked as 
well as successful jets, depending on its duration. 

We show that protons can be accelerated to 10'*-10^ 
GeV in the reverse shock. The low and high-energy pro- 
tons produce mesons efficiently by pp and pj interactions 
respectively. For energies in between, cooling mecha- 
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nisms are faster, and result in a characteristic suppres- 
sion feature in the meson spectrum. Mesons that are 
produced experience significant coohng due to the high 
proton and photon densities in the shocked jet. However, 
we stih expect a reasonable number of neutrino events. 
For instance, for a choked mildly relativistic jet with 
isotropic energy 10^^ erg at 10 Mpc, we expect ~ 2 neu- 
trino events at IceCube from the reverse shock. Interest- 
ingly, assuming the same jet Lorentz factor, we find that 
a lower isotropic energy of 10^^ erg produces more neu- 
trino events (~ 20 events). This is due to meson cooling. 
Specifically, we find that neutrinos are favored by high- 
Lorentz factor, low-luminosity, and long-duration jets, as 
meson cooling is less significant for such jets. The rate of 
core-collapse supernova within 10 Mpc is expected to be 
~l-3 yr~^ [36(1 , implying particle acceleration in reverse 
shocks can potentially be tested in the next decades. 

Neutrino emission from the reverse shock is potentially 
important for short-duration jets^ that are completely 
shocked and do not break out of the star, i.e., choked 
jets. We argue that while it is possible that internal 
shocks occur prior to jet shocking (in particular for low 
Tj jets), the allowed range of radii for internal shocks 
to occur is constrained from above by the reverse shock, 
and below by the fact that the stellar material is opaque 
to high-energy neutrinos. This suggests that neutrino 
emission from choked jets are of similar magnitude to 
our estimates. 

This paper is organized as follows. In Sec. II wc discuss 
the evolution of the jet as it propagates through the star. 
In Sec. Ill we treat proton acceleration, and suppression 
of neutrino emission due to proton and meson cooling. 
In Sec. IV we utilize the results of the previous sections 
to investigate the detectability of high-energy neutrinos 
from reverse shocks. Finally, we finish with discussions 
in Sec. V and conclusions in Sec. VI. 



II. JET MODELS 



A. Jet dynamics 



We first discuss the evolution of a relativistic Tj ^ 1 
jet as it strikes the stellar matter. The aim is to ob- 
tain the evolution of the jet head and reverse shock 
velocities. There are several numerical simulations on 
the propaga tion of an initially hot jet in the litera- 
ture [33^ lio, Ell, \^ 43 1 , as well as analytic treatments 
of initially cold jets [H 0, SI, Sll . We follow and ex- 
tend the latter. Generally, two shocks form when a jet 
strikes external matter: a reverse shock that decelerates 



^ We note that the word "short-duration" here does not mean an- 
other class of GRBs, which is often referred to as short-duration 
and hard-spectrum GRB, and is beheved to be associated with 
compact-star mergers [371. 13811 . 



the head of the jet to a Lorentz factor F/j, and a for- 
ward shock that accelerates the external material to F^,. 
The jet and stellar material moving at Th are separated 
by a contact discontinuity, and are in pressure balance. 
Equations governing the evolution of the shocks are (for 

> 1) [ia,!!! 

es/usmpC^ = F^ 1, ns/riext = 4F/i + 3, 
eh/nhTUpC^ =fh-l, n/i/rij = 4f h + 3, (1) 

where mp is the proton mass and c is the speed of light. 
Thermodynamic quantities rii and (particle number 
density and internal energy) are measured in the flu- 
ids' rest frames, and we have used the following nota- 
tions: j (jet, unshocked), h (jet head, shocked), s (stellar, 
shocked), and ext (stellar, unshocked). Lorentz factors 
are measured in the lab frame, except for Th which is the 
Lorentz factor of the jet head measured in the frame of 
the unshocked jet, 



r,r„(i 



(2) 



where [3i is the velocity as a fraction of c. The jet particle 
density is given by (for a jet with constant opening angle) 



nj(r) 



^iso 



(3) 



where Lisa is the isotropic jet luminosity. We obtain F^ 
by equating the pressure of the shocked jet head ph ~ 
e/i/3 and the pressure of the shocked stellar matter ps = 
es/3 [l6|]. From Eq. ([T]), this requires us to solve 



Tij ^ (4Ffe + 3)(Ffe-l) 
next (4f,, + 3)(f,, -1)' 



(4) 



by substitution of Eqs. ([2]) and ([3]). In what follows we 
use next from the progenitor model E20 from Ref. [io! ]. 
an initially 2OM0 rotating star that is reduced to IIM0 
by mass-loss. It has a 1.7Mq Fe core of radius rpc ^ 
2 X 10® cm, a 7.7Mq He core that extends to tho ^ 
2 X 10^^ cm, and a surrounding H envelope. However, 
the results that we present here apply qualitatively to 
other models of different masses, with or without H en- 
velopes. As an illustration of this point we include other 
progenitor models in Fig. [1] 

Due to rotation, the mass along the rotational axis be- 
comes substantially lower than the equatorial, and the 
rotational axis is an easy escape route for a jet. We in- 
ject jets at a radius r ^ 2 x 10® cm along the rotational 
axis, which corresponds to an enclosed mass of 1.6Mq. 
Since the explosion and jet mechanisms are still under 
debate, they are treated parametrically with the follow- 
ing parameters: the isotropic equivalent energy Eiso and 
duration Tj (which yield the isotropic equivalent lumi- 
nosity through Liso = Eiso/Tj), Lorentz factor Tj, and 
variability t^. We consider two types of jets, (i) mildly 
relativistic jets with F^ = 10, and (ii) ultrarelativistic 
jets with Tj = 100. 
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FIG. 1: Jet head velocity fih^h as a function of radius. 
E15 and E20 are two progenitor models of mass VoMq and 
2OM0 respectively, as detailed in the text. L51-L53 denotes 
the isotropic equivalent luminosity of the jet, Liso = 10^^- 
lO'^^ergs-V For the L53 case, two jet Lorentz factors are 
plotted: mildly relativistic (Fj = 10, solid) and ultrarelativis- 
tic (Fj = 100, dashed). For L52 and L51, mildly relativistic 
jets are shown. 



In Fig. [T] we show the jet head velocity multiplied 
by Vh as a function of radius, obtained from solving 
Eq. Q. For L\^q = 10^^ erg s""'^, we show our two jet 
types Fj = 10 and Fj = 100, for two progenitor models. 
The first progenitor model is E20 as described above, and 
the second is E15, a similarly rotating star but with ini- 
tial mass I5M0 [4^. We see that the evolution of is 
independent of Fj, and that the initial jet head velocity 
is always subrelativistic. The latter is because the ex- 
ternal density is very much larger than the jet density, 
next 3> for conceivable set of jet parameters. As the 
jet continues to propagate, the stellar density falls faster 
with radius iriext oc r~^) than the jet density (jij (x 
and the jet head effectively accelerates, reaching mildly 
relativistic values at r > tho ~ a few xlO^^cm. Al- 
though the initial value of Ph scales somewhat with Ljso, 
the radius at which the velocity reaches c is only weakly 
dependent on jet parameters. This is because the jet 
head accelerates to relativistic velocities when the exter- 
nal stellar density suddenly drops, at the edge of the He 
core. In fact, rewriting Eq. ([¥]) in the relativistic limit 
(F„F;, » 1) yields F, cx LU'r-'/^n;,'/' M- 

Note that the kinks in Fig. [T] corresponding to com- 
position changes in the progenitor model are model de- 
pendent, e.g., an abundance of oxygen resides at r 
4 X 10^ cm. However, the general features of Ph that 



TABLE I: Summary of choked jet characteristic times. Shown 
are the core crossing time tuc, and the jet duration time Tx 
required for reverse shock crossing to occur at rx = 5 x 10^" 
cm. These values apply to both mildly relativistic and ultra- 
relativistic jets. 



Liso [erg/s] 


tUc [s] 


[s] 


Eiso [erg] 


Model 




17 


12 


1 X 10^^ 


A 


1q52 


43 


34 


3 X 10^^ 


B 


10^^ 


120 


100 


1 X 10^^ 


C 



concern us, including its initially ^ 1 value and its 
reaching ~ 1 at r ~ rjjo, are model independent. We 
therefore consider only the E20 progenitor model in the 
rest of our paper. 

From the velocity profile we now obtain an estimate of 
the time taken for the jet to propagate through the stellar 
core, iHc ~ 10-100 seconds for Liso = 10^^-10^^ ergs~^ 
(see TableU]). We adopt model B in Table[T]as our fiducial 
supernova jet model. 



B. Choked and successful jets 

Since jets reach relativistic velocities as they leave the 
He core, the time taken for a jet to cross the H envelope 
tn ^ r}i/{2cTfJ <^ inc 145|], and the time taken to leave 
the star is approximately the time taken to cross the 
core, iHc ~ 10-100 seconds. If the jet duration Tj is less 
than tHc, the cold jet will be unable to break out of the 
star, i.e., a choked jet that is dark in gamma rays. On 
the other hand, if Tj > fuc, the cold jet can break out, 
and upon entering the H envelope or interstellar medium, 
receive a large velocity boost. A GRB is then possible. 

We first discuss choked jets that do not break out of 
the star. Such jets are progressively shocked and decel- 
erated by the reverse shock, until at a crossing radius 
we define as Tx, the entire jet is shocked. The reverse 
shock is relativistic in the jet frame (F^i ^ 1), and ef- 
ficiently converts the jet bulk kinetic energy to internal 
energy. Energy conversion is thus complete in just one 
jet crossing by the reverse shock P[ f48|. 

Internal shocks between jet shells of different luminos- 
ity and Tj will occur at radii no larger than Vg ~ 2r^ cty. 
If they occur before the jets are decelerated significantly 
by the reverse shock, a fraction of the jet's bulk kinetic 
energy is converted at internal shocks [50l,'5l'|. Therefore, 
the relative sizes of rs and Tx indicate at which shock (in- 
ternal or reverse) the jet's kinetic energy is converted. In 
the current work we focus on cases where < rg, i.e., 
the jet's kinetic energy is converted by the reverse shock. 
We discuss in Sec. V internal shocks occurring before re- 
verse shock crossing. 

To estimate rx for a choked jet, we note that F/i w Fj 
during propagation in the core. The reverse shock thus 
crosses the entire length of the jet in a time A/c ^ Tj, 
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TABLE II: Definitions for jets and their observables. Suc- 
cessful highly relativistic jets are the leading model for GRBs. 
The successful mildly relativistic jet and the choked jets are 
dark in gamma rays, because shocks typically occur while the 
jet is optically thick. Entry "this work" represents the specific 
focus of this work. In addition, since the reverse shock occurs 
in successful jets too, the present work applies to the entire 
neutrino (i^) row. 



Choked 


Successful 


Tj = 10 = 100 


=10 = 100 


7 rays dark dark 


typically dark bright 


u this work this work 


[32,^3] [16,64] 



where A is the jet length in the lab frame. We obtain Tx 
as the distance the jet head has traveled in time Tj . As an 
example, adopting a jet luminosity Lisa — lO^-^ergs"^, 
we obtain ^hc ~ 43 s; a jet with duration Tj — 30 s is 
therefore choked, and we find that rx ~ 4 x 10^" cm for 
both the mildly relativistic as well as ultrarelativistic jet 
cases. For calculation purposes we set rx = 5 x 10^" cm. 
This requires a jet duration of Tx ^ 0.7-0.8iHG, depend- 
ing on the value of Liso- Wc summarize duration times 
in Table m 

We note that we could have in principle chosen a 
smaller rx due to smaller Tx . However, as we discuss 
in Sec. V, the stellar matter becomes opaque to high- 
energy neutrinos for radii r < 10^*^ cm, and small radii are 
uninteresting. We also note here that rotational effects 
decrease the matter density along the rotational axis, in- 
creasing the jet head velocity at small radii and thereby 
shortening the core crossing time tne- A decrease in den- 
sity also decreases neutrino opacity. However, the effect 
is not expected to be large, because density modifications 
typically occur at radii smaller than or equal to rpe [3] ■ 

Next we discuss jets that successfully break out of the 
star, with Tj > <hc- These jets are progressively shocked 
as they move through the progenitor, but the reverse 
shock does not cross the entire jet. Therefore, for suc- 
cessful jets, only the front portion of the jet is shocked. 
We assume rx just inside of the He core radius. We 
stress that for successful jets, rx is not the point of re- 
verse shock crossing, but the last point of shocking by 
the reverse shock. 

A few words about the jet models must be noted. First, 
we have assumed a conical jet geometry in Eq. ([3]). This 
is a good approximation for jets propagating in a suffi- 
ciently sparse medium, such as the outer regions of stars. 
From analytic studies, the jet opening angle 6 oc r~^/^ for 
r < rpc and 9 oc constant for r > rpc !45| . We have there- 
fore injected jets at a radius ~ rpe with profile Eq. 
On the other hand, under conditions such as propagation 
in the dense iron core and/or when strong magnetic fields 
are associated with the launch of the jet, the high exter- 
nal matter density and/or magnetic hoop stress work to 
coUimate the jet. Numerical studies find that jets remain 
tightly coUimated in the iron core [52i] . and further be- 



yond the stronger the magnetic field [s^. For jets that 
remain strongly coUimated, the reverse shock becomes 
increasingly less important due to the more penetrating 
nature of the jet. Further, it is less likely for these jets 
to form choked jets. For illustration, consider an ex- 
treme case where a jet of luminosity 10^^ erg s~^ remains 
tightly cylindrical all the way until break out, with con- 
stant cylindrical radius 10^ cm [sl]. For such a jet, we 
find that tfjc ~ 6 s and Tx ^ 2 s. 

Second, we have assumed a type lb progenitor. A jet 
with iiso = 10^^ erg s~^ propagating through a type Ic 
progenitor would clear the core in tcNO ~ 20 seconds, 
and jets are more likely to break out. 

Third, we have defined successful jets as those which 
break out of the star. However, these do not necessarily 
always produce observable GRBs — for small Tj and t^, 
internal shocks will occur while the jet is optically thick 
to photons. For example, r^ ~ IG^^F^ i^v,-! cm (we de- 
fine Qa — Q/10" for a quantity Q in cgs units), while 
the jet typically becomes optically thin to photons at 
[(Tiso(TT)/(4^r,mpc2)]i/2 ^ IQi^^/^^^r^'y cm, where 
(Tt is the Thomson scattering cross section. The jet is 
still transparent to high-energy neutrinos. We summa- 
rize our various jet definitions in Table HIl 

To conclude, we use the following crossing radii: (i) 
rx = 5 X 10^° cm for choked jets and (ii) = 10^^ cm 
for successful jets. 

III. NEUTRINO PRODUCTION 

In this section we discuss proton acceleration and 
neutrino production through decays of pion and kaon 
mesons. The reverse shock is likely to be coUisionless, 
as we discuss in section IIIA. If first-order Fermi acceler- 
ation is realized, a fraction of the jet protons are accel- 
erated, with a spectrum 

dfip ^ Up 

where the power- law index p « 2.3 in the ultrarelativistic 
shock limit [13, H^, andp = 2 for a nonrelativistic shock. 
In this paper, we optimistically adopt p = 2 for all of our 
figures, but we present numerical results for p = 2.1, 
2.3, and 2.5 also. The minimum proton energy would 
be ep,rnin ~ T jTUpC^ , but otherwise it is not well known; 
we assume ep.min — 10 GeV, but its exact value does not 
affect our results. 



A. Maximum proton energy 

The maximum proton acceleration energy is deter- 
mined by comparing the proton's acceleration time scale 
to its energy loss time scales. Ando and Beacom have 
shown that the photopion process is the most compet- 
itive cooling mechanism in the proton energy range of 
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interest [33|. We therefore focus on photopion produc- 
tion, and refer the reader to Appendix A for treatments 
of other coohng mechanisms such as synchrotron cool- 
ing. We will however for completeness display all cooling 
mechanisms in Fig. [2] Note that the Bethe-Heitler pro- 
cess turns out to be an important cooling process for pp 
interactions, as we will discuss in the next subsection. 

The energy density in the shocked jet head is ej^ = 
i4fh+3){fh-l)njmpc' « (f ^/r,)2Li,J(7rr2c). We have 
specifically used primes (/) to remind that these quanti- 
ties are in the jet head comoving frame. We assume that 
a fraction Se — £b = 0.1 of the shocked plasma internal 
energy is converted to relativistic electrons and magnetic 
fields, in analogy to GRBs. This gives the comoving mag- 
netic field strength as 



B' « [(8eBLiso)/(r2c)]i/2(f^/rj) 



''x.10.7 




(6) 
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jet scenarios studied here. One can obtain from Eq. ^ 
that for example, F/j = 2 yields f /j/Fj k, 0.3. The rela- 
tivistic electrons immediately lose energy by synchrotron 
radiation, and the photons thermalize due to the large 
opacity n'^aT^h ~ 10^, where = (4rj -I- 3)nj(rx) « 
3 X 10^^ cm~'^ and A/j is the proper thickness of the 
shocked jet plasma. Here we adopt A/j — 0.2dr, where 
9 = 10^^ [l^. The black-body temperature is then 



T' 



[(15eeii,o?i'c3)/(^V2c)]i/4(f„/F,)i/2 



keV, 




(7) 



and the average photon energy is e'^ w 2.7T/. The aver- 
age density is 

n; « 0.33[(eeLiso)/(7rV2c2?i)]3/4(f„/F,)3/2 

/ 3/4 3/4 X - 3/2 



3/2 
''x,10.7 



The total proton cooling time is determined by the 
inverse sum t'~^l,^ = Ei t^y/ , where i runs through the 
various cooling mechanisms. As discussed in the open- 
ing paragraph, we focus on cooling by photopion pro- 
duction. Protons interact with the dense photon field 
Eq. ([5]), and for sufficient proton- photon centre of inertia 
energies mesons are produced. We determine the cooling 
time scale by t'^^ — e'p/{capjn'^A€p). We assume the con- 
ventional inelasticity K = Ae'p/e'p = [1 — (m^ — rn^)/s]/2 
where s is the invariant mass of the system. We fit (jpy 
to [5^ , and use a blackbody photon spectrum of temper- 
ature Eq. d?]). The time scale is plotted in Fig.[2l 

The acceleration time scale in the jet head frame is 
given by t'^^ = 6pR'^/c, where 9f is a constant if we 
assume the diffusion coefficient is proportional to the 



FIG. 2: Inverse of proton cooling and acceleration time scales, 
in the shocked jet head frame, as functions of the proton en- 
ergy. The case of reverse shock crossing at rx = 5 x 10^° cm 
is plotted for supernova jet model B, with other parameters 
Tj = 10, Ee = £b ~ 0.1, and = 10~^ s. Cooling mech- 
anisms shown are: synchrotron (blue, thin dashed), inverse- 
Compton (green, dot-dashed), Bethe-Heitler (green, dotted), 
proton-proton (black, long dashed), and photopion (black, 
thick dashed) . Also labeled are threshold energies f^f^' and 
^'p'th (^®^ S^*^- IIIB). Photopion limits the proton energy to 
7 X 10* GeV. 



Bohm diffusion coefficient, and i?^ is the Larmor radius. 
0p > 10 is a fairly conservative value, and 9p ^ 1 can be 
achieved for mildly relativistic shocks 01 ■ We use 9f ~ 10 
which gives agreement with previous studies [3^ . [33| . Us- 
ing the magnetic field in Eq. ([6]), we obtain 



4c « Ope'^/ieE'c 



1 X 10" 




(9) 



where is measured in GeV. 

In Fig. [21 we plot the proton acceleration and cool- 
ing time scales as functions of the proton energy, both 
in the jet head comoving frame. Acceleration is stronger 
than cooling at low proton energies, but photopion cool- 
ing becomes stronger at higher energies. It is evident 
that photopion cooling is the strongest cooling mecha- 
nism. Equating t'^^^ 
proton energy as , 

dependency is rl^lo.r£B%£e!^-iL~o!52i^Ki/Tj,i)-'^/^ 
the flat region of t'p^ . 

We note that the reverse shock is likely to be collision- 
less, because the proton plasma frequency usp ^ 10^"* Hz 



t'p p-y I we obtain the maximum 
~ 7 X ID"* GeV. The parameter 



m 



6 



is larger than the typical radiation and particle collision 
frequencies. For example, the frequency of proton-photon 
collision is Wcoii ^ n'^cTp^c ~ 10^" Hz, while the frequency 
of proton-proton collision is ~ n' UppC 10® Hz. 



B. Proton cooling and meson spectrum 

Depending on its energy, protons produce mesons via 
pp and/or p7 interactions. We consider the pion and kaon 
mesons. The multiplicity in each pp and P7 interaction is 
taken to be 1 for pions and 0.1 for kaons; this matches the 
required ratio in the relevant energies [53, Issl ls^ . Gener- 
ally, kaons produce higher energy neutrinos [33, IHSl , but 
we additionally find that they can dominate those from 
pions because of less efficient kaon cooling. 

For low energy protons, pp interactions dominate. Me- 
son production by pp interaction is efficient since the 
opacity is very high, n'^appl^h 10^, where n'^ = {ATj + 



3)nj(r> 



3 X 10^^ cm ^ and tTo„ = 5 x 10 



ing mesons since the opacity is n'^(J. 
we use (T™ — 10~^^cm^ 'S' 



^^p^L^n - 10"^, where 
For model B in Table I, 



i{BH) 



100 GeV and e 



p.th 



2800 GeV. 



We assume mesons are produced with 20% of the par- 
ent proton energy. In the absence of competing cooling 
mechanisms, the meson spectrum follows the initial pro- 
ton spectrum. On the other hand, when other cooling 
processes are faster than meson-yielding pp and pj inter- 
actions, the meson spectrum is suppressed by an amount 
equal to the ratio of time scales. We find that the Bethe- 
Heitler process is a competitive cooling mechanism for 
pp interactions, and the resulting meson spectrum will 
be suppressed by a factor given as 




r l(BH) 

otherwise. 



< el, < e: 



'(P7) 
p,th 



(10) 



We here approximate i'^^ as a power-law 
t' (e' le'^^"^)-^ for e'^^^^ < e' < e'^^^^ 



''p,th ^ . 

constant. We find that the value of k is fairly model in- 
dependent, with values k k. 0.8 (1.2) for Tj = 10 (100) 
for model B. 



C. Meson cooling and neutrino spectrum 

In this section we discuss meson cooling and the neu- 
trino spectrum. Mesons cool similarly to protons, by 
radiative and hadronic (collisions with protons) pro- 
cesses. The meson cooling time scales are analogous to 



those for the proton and can be summarized as t[^^ = 
{im^c^)/[AaTmle'{U'^ + U'^)] and t^d = e'/(ca„n;Ae'), 
where m and e' are the meson mass and energy, and 
U'^ = (-'^n'^ and U'^ — B''^/{8tt) are the photon and mag- 
netic field energy densities (we assume inverse- Compton 
process is in the Thompson regime; synchrotron domi- 
nates at high energies). Here, ah — 5 x 10^^® cm^ 56] is 
the cross section for meson-proton collisions and Ae' = 
0.8e' [l^l is the energy lost by the meson per collision. 
The total cooling time is t'~^ 
scales are 



*rad +^had' ^ud the time 



had 



3 X 10 



-7 



and 



7r,rad 



1 X 10" 



,io.7-'^i,i 



Jiso.52 



„2 

' x,10.7 



(11) 



(£e,-l + £B,-l)-^iso,52 o 



Proton cooling by pair production (Bethe-Heitler, here- 
after BH, see Appendix for treatment) overtakes pp in- 
teractions at higher energies (see Fig. [5]). We define this 
transition as At even higher energies P7 domi- 

nates over BH cooling, and we define this transition en- 
ergy as ^'p^^h ■ The P7 interaction is efhcient at produc- ^±^^± 



(12) 



for pions. For kaons, the radiative cooling time scale is 
longer, because t[^^ cx m'*. Hadronic cooling dominates 
at low proton energies, and radiative cooling becomes 
important at higher energies. 

Charged pions and kaons decay into neutrinos through 



H + Vf^{Df^), with the muon neutrino taking 
1 /4 of the meson energy. We neglect secondary neutrinos 
from muon decays since muons immediately lose energy 
by radiative cooling. We define the break energy e^^^^j^ by 
equating j't — t'^ ^ ^had' where 7' and r are the meson 
Lorentz factor and proper lifetime. Below e'-^^^ , neutri- 
nos have a typical energy 0.056^ and a flat power per 

decade neutrino spectrum. Above e'^l^, the spectrum is 

suppressed by a factor thad/(7'''')- Next, we define e^^^^ 
when radiative cooling begins to dominate, by equating 
*had = *rad- ^'^^ ^o rclativistic cffects, the neutrino en- 
ergy in the observer frame is related to the parent me- 
son energy in the jet head frame as ~ rfte5j/4 and 
e„ — Thex/"^- We define the suppression function C as 



1 



for < e 



C(e.) 



(1) 

brk/ 
(1) ,(2) / 2 



€1 — \l 



for e 



(1) 



(1) 



for > t^^l^^. 



— '^^^ ''i/.brk 



(13) 



For some supernova jet parameters, it is possible that 
meson goes from decay dominated straight to radiation 
cooling dominated. For these cases we define the break 
energy e^'^J^ by equating 7'r = and the suppression 
function C as 



C(e.) 



1 



for < e 
for £,, > e 



(r) 



(r) 
■,brk- 



(14) 



The break energies for neutrinos from pion decay, in a 
choked mildly relativistic jet, are 

,7^(2) 



<brk = 0.4GeV 



= 10 GeV. 



(15) 
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The parameter dependencies are io.7rj,ir/i,o-^iso 52 
and r^,o(ee,-i + £i3,-i)"^r-i(f /..i/r^- 1)-^ for the first 
and second break energies, respectively. On the other 
hand, break energies for neutrinos from kaon decay are 
higher due to the higher kaon mass and faster decay time; 



K{1) 
f,brk 



5GeV 



K(2) 



3000 GeV. 



(16) 



The parameter dependencies are the same as for pions. 
The break energies are smaller than those of internal 
shocks [33} with the same jet parameters, by as much 
as ~2 orders of magnitude. This is a combination of the 
higher proton and photon densities in the shocked jet 
causing more effective cooling, and the fact that the neu- 
trino energies in the shocked jet head receive no Lorentz 
boosting. 

Next, the break energies for a choked ultrarelativistic 
jet are 



K{1) 
-i/.bik 



60 GeV 



'=iy,brk ~ 
K(2) ^ 
i/,brk 



0.9 GeV, 
= 300 GeV, 



(17) 
(18) 



which can also be obtained from scaling Eqs. P^ - P^ by 
their dependencies on T j Importantly for neutrino emis- 
sion, suppression due to radiative cooling is independent 



(1) 



of Tj, since the product of 

,brk iy,brk' 



brk 



X e 



(2) 
.brk 



is unchanged, 
indicating that radiative 



Also, note that e 
cooling dominates the cooling process for pions. 
case, 



In this 



ti/,brk ^ ^ 



,10 



SB. 



Nl/2 7-1/2 




GeV, 
(19) 

We show in Fig. [3] the resulting suppression factors 
for neutrinos from kaons, for our mildly relativistic and 
ultrarelativistic jets. Suppression is less intense for rel- 
ativistic jets. We also shade in for reference the energy 
range where the meson spectrum is suppressed by Bethe- 
Heitler cooling, from 5-140 GeV. Note that for pions, 
radiative cooling dominates for energies of interest. 



IV. DETECTION OF NEUTRINO BURSTS 

In this section we estimate the emission and detectabil- 
ity of high-energy neutrinos from reverse shock acceler- 
ated protons. We consider a supernovae at a distance 
oi Dl, and use the code ANIS (All Neutrino Interaction 
Generator) [b^I to calculate the neutrino induced muon 
spectrum at a neutrino detector with an effective area 
of 1 km^, i.e. an IceCube class detector utilizing upgoing 
muons. First, we estimate the fluence of neutrinos as ^33] 



p,min) 



(20) 



where (n) is the meson multiplicity (1 for pions and 
0.1 for kaons), is the branching ratio of meson de- 
cay into neutrinos (1 for pions and 0.6 for kaons), and 




100 1000 
Neutrino Energy, e^^ [GeV] 



10" 



FIG. 3: Suppression factor as a function of neutrino energy, 
for neutrinos from kaon decays. Shown are the choked mildly 
relativistic jet (solid lines) and choked ultrarelativistic jet 
(dashed lines). The kaon experiences less cooling for an ul- 
trarelativistic jet. Note that pions are always in the radiation 
cooling regime for neutrino energies shown. We also shade 
in for reference the energies where Bethe-Heitler cooling is 
competitive and suppresses the meson spectrum. 



is the fraction of proton energy carried by neutri- 
nos in the absence of energy loses; 1/8 for pions and 
1/4 for kaons, since neutral and charged mesons are pro- 
duced with roughly equal probability, and muon neutri- 
nos carry roughly 1/4 (1/2) of the pion (kaon) energy in 
meson decay. The functions CBHi^iy) and C(ci^) are the 
suppression factors due to proton and meson coolings, 
respectively, and the ln(ep nj^x/^p min) factor normalizes 
the proton spectrum to the jet energy (for p = 2). 

In Fig. |3]we show the cumulative spectrum of neutrino 
induced muons, from a supernova with a choked jet. We 
show three cases - a supernova distance of Dl = 10 Mpc 
possessing a mildly relativistic Tj — 10 jet (solid lines), a 
supernova distance of Dl — 10 Mpc possessing an ultra- 
relativistic Tj — 100 jet (dashed lines), and a similarly ul- 
trarelativistic jet but supernova distance = 400 Mpc 
(dashed lines labeled 400 Mpc). We use model B with 
parameters £e = Sb = 0.1 and = 10~^ s. We take 
into account the muon range, which effectively increases 
the detector volume. We evaluate the muon energy 
when it enters the detector if it is produced outside, or 
at the production point if it is inside. Considering a 
detection threshold muon energy of 100 GeV, we expect 
about 6 muon events from a 10 Mpc supernova possess- 
ing a choked mildly relativistic jet. Similarly, we expect 
about 28 from an ultrarelativistic jet. We have plotted 
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FIG. 4: Cumulative muon event number from choked jets 
with model B parameters. Thin-black (thick-green) lines de- 
note pion (kaon) contributions. The majority originates from 
kaons. Solid lines denote a mildly relativistic jet and a super- 
nova distance of 10 Mpc. Dashed lines denote an ultrarela- 
tivistic jet and two supernova distances - 10 Mpc for illustra- 
tive purposes and 400 Mpc for event calculations, as labeled. 
Due to weaker kaon cooling, an ultrarelativistic jet yields more 
neutrinos than a mildly relativistic one. The spectrum be- 
comes almost flat around ~ 100 GeV, due to Bethe-Heitler 
suppression of the meson spectrum. The total number of 
events expected above the detector threshold energy is about 
6 for the mildly relativistic case. For the ultrarelativistic case 
with distance 400 Mpc we expect 0.02 events. Also shown is 
the atmospheric neutrino background over 1 day in a 3° circle. 



the ultrarelativistic jet case at an optimistic distance of 
10 Mpc to illustrate the effect of meson cooling. Although 
the pion contributions are similar between the two jets, 
the kaon contribution is larger in the ultrarelativistic case 
due to the weaker suppression. In the remainder we as- 
sume the distance 400 Mpc for ultrarelativistic jet calcu- 
lations. For model B, this yields 0.02 events. 

Note that since the muon spectrum is steep, a lower 
detector energy threshold will yield much higher events. 
Additionally, one may observe the suppression due to 
Bethe-Heitler cooling, which is expected to manifest be- 
low 0.05e'^^J^ 140 GeV. For a steeper initial pro- 
ton spectrum, motivated from studies of ultrarelativis- 
tic shocks, the expected event is 4 (2, 1) for p = 2.1 
(2.3, 2.5). In all jets, the majority of neutrinos come 
from kaons, because kaons are more massive than pions 
and suffer less cooling (the radiative cooling time scale is 
oc m^). Also, the kaon decay time is slightly shorter than 
the pion. The events discussed here arrive in a « 34 sec- 
ond time bin and ~ 3° angular bin, allowing very strong 



FIG. 5: Same as Fig. |4] but for model C parameters. Ap- 
proximately 17 (0.04) events are expected above the detector 
threshold for a single mildly relativistic (ultrarelativistic) jet 
and supernova distance 10 Mpc (400 Mpc). 



rejection of atmospheric neutrino backgrounds. We show 
in the figure the atmospheric neutrino background within 
1 day and a 3° bin for comparison. The expected neu- 
trino signal from a choked mildly relativistic jet easily 
exceeds the atmospheric neutrino background. 

Next we discuss a scan of supernova and jet param- 
eters. First, as shown above, a larger Tj reduces kaon 
suppression and increases neutrino emission. Next, we 
consider jet model A of Table HI i.e., with larger Liso but 
smaller Tj, and an overall larger Eiso- Since the neutrino 
fluence scales with Eiso, one naively expects more neu- 
trinos. However, the expected muon event rate is in fact 
smaller, 2 (0.01) for a mildly relativistic (ultrarelativistic) 
jet. This is because of stronger meson cooling: the proton 
and photon densities in the shocked jet scale as rip cc Liso 

and n'^ oc L^^^ , so a larger Liso results in stronger me- 
son cooling. We conclude that neutrino emission from 
choked jets are more favorable from low-luminosity and 
long-duration jets, even if .Eiso is smaller (assuming the 
same jet Lorentz factor). Indeed, model C predicts 17 
(0.04) events at IceCube. We show the cumulative muon 
event spectrum for jet model C in Fig. [S] 

In Fig. [6] we show the same as Fig. |4] but for successful 
jets. We show a mildly relativistic Tj = 10 jet at 10 Mpc 
(solid lines), and an ultrarelativistic Tj — 100 jet at 
400 Mpc (dashed lines). The other parameters used are 
Lisa = lO^^grgg-i^ r^. ^ iQQg^ ^ o.ls, £e = es = 0.1, 
and Tx — 10^^ cm. With these parameters, ino ^ 43 s, 
and the cold jet does indeed break out of the star. We in- 
troduce an additional factor of 1 /2 to the neutrino fluence 
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FIG. 6: Same as Fig. |4] but for jets with successful jets. 
The supernova and jet parameters are Lino = lO^^ergs"^, 
Tj — 100 s, and —0.1 s. About 36 (0.1) events are expected 
above the detector threshold muon energy from a mildly rel- 
ativistic (uhrarelativistic) jet supernova at a distance 10 Mpc 
(400 Mpc) away. 



for successful jets, because only about ^ t-nc/Tj ~ 1/2 
of the jet length is shocked by the reverse shock. Note 
that our chosen parameters have the equivalent Eiso to 
jet model A. However, we expect - (1/2) (22/10-1) ~ 20 
times more neutrino events than jet model A, because 
iiso is a factor 10 smaller and is a factor 2 larger. 
Indeed, evaluating the muon event rate yields about 36 
(0.1) events per mildly rclativistic (ultrarelativistic) jet. 
As in the case of choked jets, most of the contribution 
comes from kaons. These events will show temporal and 
angular clustering, and in addition, these will appear as 
precursor signals to GRB photons, allowing very strong 
rejection of background. 



V. DISCUSSION 

Our result is that neutrinos from the reverse shock are 
strongly suppressed due to proton and meson cooling. 
However, we still expect a reasonable number of events 
at neutrino detectors. One reason for this is that pro- 
ton cooling affects only energies near and below current 
detector threshold energies. This produces a suppres- 
sion feature in the spectrum but does not decrease event 
counts at IceCube. For our models A, B, and C in Ta- 
ble HI we expect 2, 6, and 17 muon events respectively at 
a IceCube class detector. Here, we have calculated for 
a supernova 10 Mpc away, possessing a Tj = 10 mildly 



rclativistic jet, and a p = 2 initial proton spectrum. If 
we assume a steeper index, p — 2.1 (2.3, 2.5), we expect 
4 (2, 1) events for jet model B. In all cases, the neutrino 
flux is dominated by decay of kaons, since cooling is less 
efficient for kaons than pions. 

Since neutrino emission from the reverse shock is sup- 
pressed by meson cooling, we find that neutrino emission 
is favored by parameter sets for which meson cooling is 
weaker - high-Lorentz factor, low-luminosity, and long- 
duration jets. This fact is readily observable from the 
neutrino event predictions of models A, B, and C. As fur- 
ther illustration, if we assume a Vj = 100 ultrarelativistic 
jet we expect 11, 28, and 55 events from jet models A, B, 
and C respectively. Note that supernova and jet param- 
eters are expected to show considerable scatter. Hence, 
given neutrino emission is cooling suppressed, the event 
count need not necessarily correlate with the total burst 
energy. 

Neutrino emission from the reverse shock is partic- 
ularly important for choked jets, which are completely 
shock decelerated by the reverse shock. Is it possible for 
choked jets to produce more neutrino emission than our 
above estimates? If the duration is long enough such that 
Ts < rx, it is possible that internal shocks accelerate pro- 
tons in the unshocked part of the jet, while the jet head is 
being shocked. The internal shock accelerated protons do 
not travel far, producing mesons efficiently via pp interac- 
tions with the large target jet proton density J52, 33]. For 
high-Lorentz factor jets, the protons can travel and inter- 
act with termination shock photons In both cases 
meson cooling is less severe than in the reverse shock, and 
we can expect more neutrino emission. However, internal 
shocks must occur at large enough radii that the stellar 
material is not opaque to high-energy neutrinos. The 
neutrino opacity is — Next(i')<Jvp(ejj), where N^^xt is 
the column number density of the star and aup is the sum 
of neutrino charged-current and neutral-current cross sec- 
tions, which is approximately oc in the energy range in 
question. For a 10^ GeV neutrino emitted at r = lO^" cm. 
Next « 8 X lO^s cm-2 j^j^j cr,,p = 9 X 10-37 cm^ and 
we find that t^, « 0.68. For a 10'^ GeV neutrino, « 6. 
Thus, only neutrinos with energies ti, < 10^ GeV can 
leave the star from r = 10^*^ cm. Clearly, these results 
depend on the stellar density profile assumed, and rota- 
tional effects have been neglected in our simple estimate. 
However, it seems clear that the stellar matter is opaque 
to high-energy neutrinos for small radii, constraining the 
range of allowed radii for internal shocks. 

The case for the reverse shock in successful jets, which 
break out of the star, seem better. This however is not 
due to any intrinsic property of successful jets. It is 
merely because for a given total jet energy, successful 
jets have longer durations and lower luminosities com- 
pared to choked jets. This results in less meson cooling; 
for example, we expect 2 events from a choked jet (of 
energy Eiso — 10^^ erg), while we expect 36 from a suc- 
cessful jet of the same total energy. However, it should 
be noted that we expect internal shocks to be impor- 
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FIG. 7: Suppression factor for neutrinos from the reverse 
shock (solid lines) and the internal shock (dot-dashed lines). 
We adopt our successful mildly relativisitc jet (Liso ~ 
lO^^ergs"^, Tj = 100s, ty — 0.02s), and show the suppres- 
sion factors for pions (thin-black lines) and kaons (thick-green 
lines). The vertical dotted lines indicate break energies. 



tant for successful jets. Considering our successful jet 
with Liso = 10^52 erg s-\ Tj = 100 s, and ty = 0.02 s 
(so that internal shocks occur at Ts « rx), we find that 
the total muon event number from the internal shock is 
'^lOO times larger than our estimates from the reverse 
shock [ll, m, m, [131 ■ This is due to less meson cooling 
in the cold jet, and we compare in Fig.[7|the suppression 
factor for mesons in the cold jet and hot jet. It is clear 
that suppression in the cold jet is much weaker. 

Despite being overwhelmed in terms of numbers, re- 
verse shock neutrinos will arrive prior to neutrinos from 
internal shocks. Meszaros and Waxman [l^ studied pre- 
cursor signals arising from internal shock accelerated pro- 
tons interacting with photons at r ~ 10^^ cm, while Raz- 
zaque, Meszaros, and Waxman [g^l extended this to in- 
clude r ~ ruo. Neutrinos from reverse shocks will pre- 
cede even these precursor signals, provided the radius is 
large enough that the stellar material is thin to neutri- 
nos. The temporal profile of neutrinos therefore offers 
the possibility to distinguish locations of particle accel- 
eration for exceptionally close supernovae. The expected 
time difference between components is of order a few to 
tens seconds, which is within the time resolution of neu- 
trino detectors. 

We classified jets by whether they break out of the 
star or not, depending on their durations. Choked jets 
are bright in neutrinos and dark in gamma rays, while 
successful jets are bright in neutrinos but may or may not 



be bright in gamma rays, depending on whether internal 
shocks are optically thick or not. A successful yet gamma 
dark jet may be typical for mildly relativistic jets, where 
the smaller Tj implies internal shocks typically occur at 
smaller radii (32, . . These could be distinguished from 
choked jets by the presence of an afterglow like radio 
signature, for example. 

In the current work we have assumed a type lb progen- 
itor. However, if the progenitor is a type Ic, the situation 
is slightly different. Having no He core, the core crossing 
time is shorter, icNO ~ 20 s, for Ljso = 10^^ erg s~-^, 
and choked jets will become rarer. This implies the rel- 
ative importance of internal shocks over reverse shocks. 
For successful ultrarelativistic jets in type Ic progenitors, 
the neutrino precursor signal can be as strong as those 
discussed previously ^16i , j64| . Our results will however 
apply well to type II progenitors. 

The core-collapse supernova rate within 10 Mpc is '^1- 
3 yr~^ [3^], and therefore, one could test the possibihty 
of particle acceleration in the reverse shocks in the next 
decade. In particular, correlating with the optical signal 
of the supernova explosions (not necessarily the jet-like 
signature) [i^ or MeV neutrino detections [3^ would 
greatly help enhance the significance of the signal. Ex- 
amples of recent close supernovae include SN 2007gr in 
NGC 1058 (distance -10 Mpc) [11], and a further Hst 
can be found in Ref. [H^]. The actual rate could in fact 
be higher, as a systematic survey of all nearby galaxies 
has not been performed to date. 



VI. CONCLUSIONS 

We showed that neutrino emission from reverse shock 
accelerated protons are strongly suppressed due to meson 
cooling in the hot jet and hot stellar material surround- 
ing the reverse shock. Despite this suppression, we still 
expect a fairly promising number of events at IceCube, 
ranging from 2 to 20 for a close supernova of 10 Mpc host- 
ing a choked mildly relativistic jet of total isotropic en- 
ergy 10^^-10^'* erg. We find that for a supernova of given 
energy i^iso , neutrino emission is favored by high-Lorentz 
factor, low-luminosity, and long-duration jets, for which 
meson cooling is less severe. 

Our result sets the scale for neutrino emission from 
choked jets. For choked jets to provide more neutrinos, 
internal shocks are required to occur in the unshocked re- 
gion of the jet. However, this is geometrically constrained 
from above by the reverse shock, and from below by the 
fact that the stellar matter is opaque to high-energy neu- 
trinos for small radii. 

Regarding neutrino precursor signals which precede 
GRB photons, we conclude that the reverse shock neutri- 
nos likely plays a small role. However, they will precede 
those neutrinos from internal shocks, allowing the tempo- 
ral profile of precursor neutrinos from a close supernova 
to provide an identification of proton acceleration at the 
reverse shock. 
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APPENDIX A: PARTICLE COOLING 
MECHANISMS 

We discuss cooling mechanisms in the shocked jet 
head, which act on the proton (Section IIIA, shown in 
Fig. [2|) and mesons (Section IIIC). We exphcitly denote 
by primes (/) the jet head comoving frame. All numerical 
values are derived for the proton. 

First, the same magnetic field responsible for electron 
synchrotron radiation also causes energy losses in protons 
and mesons, given by 

t'sync = (67rm;^C^)/(0-TTOeB'^4) 

- 0.04 ( 'Y"' , ] f^)'s, (Al) 

where me is the electron mass, and tux and are the 
particle mass and energy. The numerical value is derived 
for the proton. Due to the strong dependency on the 
mass, oc m^, synchrotron cooling is most effective for the 
pion meson and least effective for the proton. 

Second, protons and mesons lose energy by inverse 
Compton (IC) scattering. The most natural target pho- 
tons are due to radiation cooling of ultrarelativistic elec- 
trons. Although a full treatment using a black-body 
spectrum is simple, a simplistic approximation using 
Eq. ([7]) and Eq. ([5]) is sufficient to demonstrate that 
IC cooling is subdominant. Below a transition energy 
nipC^/T!^ ^ 7 X 10"* GeV (for the proton), inverse Comp- 
ton scattering proceeds in the Thomson regime, and the 
cooling time scale is given 

t'icTh = (3m;^c^)/(4crTTO^e^r^n^) 

- 0.04 { , \ (^)'s.(A2) 



Again, note the strong dependency on the particle mass. 
Above the transition energy, inverse compton scattering 
proceeds in the Klein-Nishina limit, where the cooling 
time scale becomes 



Third, we consider cooling by the Bethe-Heitler (BH) 
interaction (p7 — > pe"'"e~), which is important for pro- 
tons. Our treatment is the same as that described in 
Ref. [s^l. The BH cross section rises logarithmically 
with energy, obh ar\{(2%l^)\T{{2f'pf'^)l{mj,m^c^) - 
106/9]} [111- The produced are at rest in the cen- 
tre of mass frame of the proton and photon, and so the 
energy lost by the proton in each interaction is given 
by Ae^ = IvrieC^i^ ^ , where 7^ „ = (e^ + e'^)/{mlc^ + 
ie'^e'^Y/"^ is the Lorentz factor of the centre of inertia 
as seen in the comoving frame. The energy loss rate 
is de'pjdt' = n'^cGBH^^'p, and thus the cooling time is 
^'bh = ^'p/{de'p/dt') = e'p/ {2n'^caBHmeC^y^^). 

Lastly we discuss cooling due to collisions with pro- 
tons. Assuming 20% of the proton energy is lost in each 
collisions with a proton, the cooling time scale is 



tpp = <^'p/{c(Jppn'p/\e'p) 

« 1 X 10-*^ { "IimEiA s, (A4) 

\ -^iso,52 J 



for the proton, where Gpp = 5 x 10^^® cm^ [sHi. Anal- 
ogous equations hold for mesons, e.g. Eq. (jlip for the 
pion. 
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